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^- Abstract 



Collective interactions of a beam of neutrinos/antineutrinos traversing a dense mag- 
netized plasma of electrons/positrons, protons and neutrons are studied with par- 
S^ . ticular reference to the case of a Supernova. We find that the ponderomotive force 

H \ exerted by neutrinos gives, contrary to expectations, a negligible contribution to 

the revival of the shock for a successful Supernova explosion, although new types of 
convection and plasma cooling processes induced by the ponderomotive force could 
be, in principle, relevant for the dynamics itself. 
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At the end of their hves, massive stars (M > 8Mq) develop into Supernovae of type II 
and, after explosion, neutron stars are born. The mechanism of type II Supernovae is 
approximately well understood, although some obscure points remain |l|, p|. 

The iron core of the progenitor star has a mass around the Chandrasekhar limit 
(~ IAMq) so it has no stable configuration and will collapse. This is also accelerated by 
the negative nuclear pressure. During the collapse, the key process is the electron capture 
by free protons and transitions in complex nuclei, producing (electron) neutrinos that 
escape from the star, thus deleptonizing it. Electron capture would continue indefinitely, 
and the final electron fraction Y/ at collapse would be extremely small, were it not for 
the trapping of neutrinos taking place when the density in the core has reached very high 
values (p ~ 10^^ -^ 10^^ g/crn?). After neutrinos are trapped, due to (weak neutral current 
mediated) elastic scattering by nuclei, they fill available phase space and their distribution 
can be approximately described by a Fermi-Dirac distribution with effective temperature 
Tjj and chemical potential ^y (different from those of electrons). The collapse can be 
stopped only when the nuclear pressure becomes positive, which happens when the nuclei 
touch and fuse together, forming nuclear matter. Thus, collapse continues until the central 
density becomes substantially greater than nuclear density (~ (3 ^ A)^l{]^^ g / cm?) . At 
this point, the inner core rebounds and sends a shock wave out into the star. The shock 
forms not in the centre of the star, but near the surface of the homologous core (slightly 
outside it). 

In few cases, if both the temperature at the beginning of the collapse and the mass of 
the iron core are low |0|, the shock will go through the entire star and expel most of it, 
thus giving a Supernova {prompt mechanism) . The core remains behind and will become 
a neutron star; the large negative gravitational energy of that neutron star provides the 
energy to expel the mantle and the envelope against gravitational attraction of the core. 
The major part of the released gravitational energy of the neutron star goes into the 
emission of neutrinos, which mainly occurs after the Supernova material has been set in 
outward motion. 

Behind the shock, nuclei are dissociated into nucleons at a cost in energy of about 8.8 MeV 
per nucleon, and this gradually drains the energy of the shock. If the amount of material 



the shock has to traverse from the homologous core before it emerges from the iron 
core is not small, the prompt mechanism fails, because the shock rapidly looses energy 
to dissociate nuclei into nucleons. The shock is therefore apt to stall at some point 
(typically at r ~ 400 Km) and then it turns into an accretion shock in which additional 
infalling matter accretes to the existing core. The outward motion has then stopped 
and the prompt shock has failed to expel the outer part of the star. Furthermore, once 
the shock slows down, there is time for copious neutrino emission, which further saps its 
energy. However, if the shock turns into an accretion shock, neutrinos from the core can 
be absorbed (after about 0.5 4- 1 s) by material at r ~ 100 -^ 200X171, and can heat this 
sufficiently to revive the shock, which will then expel the material from the star {delayed 
mechanism) |Q. 

However, in most cases, the outgoing material (ejecta) is too far away to effectively absorb 
the neutrinos from the core but, in numerical simulations, it is found that between ejecta 
and the core there is a low density region (p ~ 10'^ g/cm^) filled mainly by radiation 
(the radiation bubble) [|]. In this bubble, e"*" e~ pairs can absorb neutrino energy by 
elastic scattering (also uu-^ e+e^ can occur), thus it continually receives new energy and 
continually exerts pressure on the ejecta, then driving the shock. We have to note that 
density increases outward near the bubble outer surface, so that dilute material in the 
bubble has to support and push down material in the ejecta. This causes Ray leigh- Taylor 
instability and convection will occur. In the delayed mechanism, convection is very fast (it 
is practically instantaneous compared to diffusion time scales) and reduces the energy loss 
by re-emission of neutrinos. Furthermore, there is an important effect. The absorption 
of neutrino energy takes place mainly at r ~ 100 -r- 200 Km; the convection brings this 
energy to the shock front, whether that is at r ~ 300 Km or r ~ AOOO Km,, and thereby 
continually supports the shock. This will succeed in ejecting the material outside the 
bubble. 

These are the essential features of Supernova dynamics, although there are still some 
doubts especially on how the explosion mechanism can undergo effectively. In particular, 
it is clear that neutrinos are responsible for transporting energy from the central core to 
the layers surrounding it, but it seems that the coUisional loss rate is marginal to produce 



the required heating. Moreover, besides the key role played by neutrinos, a complete 
understanding of hydrodynamic instabilities and overturn before, during and after core 
collapse is also crucial. 

In this paper we will concentrate on the possible collective interactions between neu- 
trinos and the dense plasma in connection with the problem of the heating of the shocked 
envelope in a Supernova, as proposed in [§],0. The idea is that the material in the shock 
is continually run over by an intense neutrino flux, so that stimulated scattering processes 
can take place, in analogy with stimulated Raman or Compton scattering for laser cou- 
pling to plasma oscillations. Thus we first study the (weak interaction) ponderomotive 
force exerted by neutrinos on the background electrons, positrons, protons and neutrons, 
and then apply the results to the case of Supernova explosion. 

Let us consider neutrinos with momentum p^ whose phase space distribution is f^(t, r, p,^) 
((t, r) is the 4-position), acting on the background particles a = e^,p, n of a plasma with 
momentum p^, distributed according to /^(t, r, p^), through an effective potential 

(the factor 2 is the statistical weight for the fermions a = e^,p,n). The force exerted 
by neutrinos over (a single) background particle a = e^,p,n can then be written in the 
form IHI 



n^ J (27r)3 J (27r)3 



^u = -2 / ^ / ,^ (/.(t,r,p.)VK//(t,r,p„) + 



- V (Ut,r,p,) (^] Ut,r,p^)\) (2) 

where the particle number densities are 

n„(t,r) = 2 y — -yj^(t,r,p^) (3) 

Mt,r) = J — ^/^(t,r,p^) (4) 

(neutrinos are completely polarized particles). We now have to calculate the effective po- 
tential for u — a coherent interactions. Let us first note that in Supernovae large magnetic 



fields B are, in general, present so that neutrino interaction with the background is 
modified with respect to the case in which the medium is non magnetized. However, since 
the calculation of the effective potential proceeds from that of neutrino self-energy in the 
considered medium [|1^, ||lT|, |[1^, in general one can write 



where the first term refers to the B = case while the second one explicitly involves the 
presence of the magnetic field. 

Let us first examine the contribution independent on B 0. Two kind of Feynman diagrams 
contribute to neutrino self-energy [0: one involves charged current interaction between 
Ve and electrons while in the others we are concerned with neutral current interactions 
between each neutrino flavour and electrons, protons and neutrons. For neutrino-electron 
interactions we have 

K}/ = V2Gp(7y/e-(t,r,Pe-) (6) 

where Gp is the Fermi coupling constant and 
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(9) 

with 6w the Weinberg angle of the electro-weak standard model. The effective potentials 
for neutrino-positron and neutrino-proton are instead given by 

V:^f = -V2GFgvfeAt,r,P.+ ) (10) 

V!ff = -v^Gp^^^/,(t,r,pp) . (11) 

Finally, neutrino-neutron interactions are described by 

K// = -^9vfnit,r,p^) . (12) 



^In what follows we consider only isotropic media or, more in general, anisotropic ones with equal 
fluxes of electrons moving in opposite directions (see |l2[ ). 
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Note also that the effective potential for antineutrino interactions is opposite to that for 
neutrino. Then, from Eq.(|^), we obtain the B = contribution to the ponderomotive 
force experienced by e^, p, n due to a beam of neutrinos/antineutrinos: 

"fJ + Qv'^Vin, -n^)) (13) 

(14) 
(15) 

{riu = n^, + ni>^ + ni>^). 

Let us now consider the contribution induced by the presence of the magnetic field (the 
second term in (^), which we suppose to lie in the positive z direction. For neutrino- 
electron interactions, by indicating with \z the polarization of the electrons in the plasma, 
we have [0 

V^^f = - V2 Gf gA cos as Xzfe{Pez,n,\z) (16) 

where cos as = Pi^-B and 

(17) 
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(19) 

In (|T^), fe{Pez,n, Xz) is the electron distribution function (for simplicity we have sup- 
pressed the spatial dependence) which for a magnetized medium takes the form 

/.fe.,«,A.) = (l + exp{ ^'^'""-^^-'-^ ])"' (20) 

where e{pz, n, A^) are the quantized Landau energy levels given by (e > 0) 



e(pe.,n,A,) = ^pI + ml +eB{2n+ 1 + A,) (21) 



with n = 0, 1,2, ..., A^ = ±1. For neutrino-positron interactions the same form in ([Tq ) 
holds, but with the replacements K//^ — Vg// and fi-^ — /i. Instead for neutrino- 
nucleon interactions the effective potential is still proportional to the polarization of the 
considered background particle, but now, for the case of interest in a Supernova, in which 



nucleons are strongly (or not weakly) non degenerate, this will be given by the proton 
fip = 2.79e/2MM or neutron /i„ = —1.91e/2MM (Mat is the nucleon mass) magnetic 



moments. We therefore have 112 



^eff = -^kT^ ^v -7F- cos aB (22) 



where C^jCy ~ 1.26 is the axial to vector nucleon form factor 



Vr// = ^^n„^cosas (23) 



Let us now pass to calculate the ponderomotive force induced by the effective potential 
in (|T6|); first observe that, since the energy levels are quantized, we have 



(2vr)3 ^^ ^ y_^ (27r)3 



tH^^^Bll^^ . (24) 



Then from Eq. (|^) we get 

Tl = V2GFgA Y. H^TT— — ^ J dpez Xzfe{Pez,n,\;,) J —^ COS Ofi V /j. 



n=0 A^ 



(25) 



Since fu does not depend on neutrino angles, 

J (T^s ^°^ (^B^^ fu = 7^-^ / vldpyV fy I dcos as cos a^ = (26) 

and so 

K = . (27) 

In a similar fashion, we can also prove that 

J^f = Tl = Tl = ^ . (28) 

We then have the important result that the magnetic field dependent term in the effective 
potential (^ does not contribute to the ponderomotive force of neutrinos by virtue of the 
dependence of Vgjj 7^ on cos as = Pi^-B, and thus Eqs. ([T3|)-(|T5D are general results. 



Some considerations are now in order. Due to weak couplings, we have that, as regards 
nucleons, 

l^^l ^ 8% I J"; I (29) 

so the neutrino force on neutrons is more efficient than that on protons; secondly, for a 
given neutrino net number, the neutral current contribution to J-"^ is about 4 % of that 



of charged current one. 

Considering the typical situation present in a Supernova, we can also predict the sign of 

the ponderomotive force exerted by neutrinos. For z/^ , Vj. the neutrino-nucleon elastic 

scattering cross section is (slightly) higher than the one for antineutrino-nucleon [p!3[| , 

so antineutrinos escape more easily from the star, leaving the Supernova more rich in 

neutrinos: 

nu,,. - n^,,. > . (30) 

The same approximately holds for the thermally produced z/g ,Ve; moreover there are also 
Ue from the deleptonization of the Supernova core, and then 

riu, - ^77, > . (31) 

Then, the net neutrino number density is a decreasing function of the distance r from the 
centre of the star, so that its gradient is negative: 

V (n,, - n^J < . (32) 

Summing up, we found that neutrino ponderomotive force is negative (i.e. attractive 
towards the centre (oc — r), assuming spherical symmetry) on neutrons and positrons, 
while is positive (i.e. repulsive from the centre) for protons and electrons. 

Such a ponderomotive force exerted by neutrinos can play, in principle, an important 
role in Supernova dynamics. The major effect would be related to the revival of the shock 
by means of the energy deposited in the background material behind the shock itself. In 
the standard delayed mechanism, if a is the averaged cross section for neutrino-electron 
and neutrino-nucleon scattering, 

(33) 
(34) 
(35) 

(where fc is a given constant depending on the particular channel, while uJi, 002 are (smooth) 
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functions of neutrino energy Ey] see Ijl^ for more details), the energy gain of a given 



particle in the background at a distance R is [|T| 

dE Ca 



(36) 



dt AttR^ 

C being the luminosity in neutrinos or antineutrinos. This relation is not, in general, 
complete since the material absorbing the neutrinos can also emit neutrinos spontaneously; 
the total energy change is then given by (PB|) times a correction factor estimated as M 



^ - O iB 

where R^- T^, are the radius and temperature of the neutrinosphere (while T is the 
material temperature). For typical numbers, 

C ^ 5x10^"^ ergs /s 
<El> ^ QTl ^ 100 MeV (38) 

R ^ 150 Km (39) 

and assuming that all nucleons are free behind the shock, the absorption average energy 
gain per nucleon is (from (|36|)) 

— ~ 50Mev/s (40) 

while for neutrino-electron scattering the rate is about 1/2 of this, due to the different 
cross section. 

However, there are also collective interactions of neutrinos on the background plasma 
described by the ponderomotive force in ([l3|)-(p!5D. The energy change induced on e^,p, n 
in the material by this force is obtained from 

where Vq, is the velocity in the medium of the considered particle. By comparison of 
(PI) and (^) we immediately see that, as expected f\, collective effects are in general 
more relevant that incoherent ones since, while the rate in (^) is proportional to Gl^, the 



expression in (O) is of order Gp and then the mechanism more efficient. However, for the 



^In |7| there is an overestimation of the momentum transfer from the ponderomotive force to the 
background, leading to exceedingly large values for the neutrino force. 



present case, there are also two important suppression factors to be taken into account, 
namely the dependence of J> on the difference between neutrinos and antineutrinos num- 
ber densities (which would be almost equal for z/, z7 produced thermally by e+e"^ uu) 
and the dependence od dE/dt on the average angle between the ponderomotive force 
and the velocity of the background particle. Let us first calculate the maximum energy 
gain implied by (0); to this end we have first to estimate the net neutrino number in 
a Supernova. We adopt a very simple model in which spherical symmetry is assumed: 
neutrinos are emitted, with a Fermi-Dirac distribution of temperature T^, from the neu- 
trinosphere (practically coincident with the core surface) located at R^ and escape freely 
outward from that. No neutrino source is considered outside the neutrinosphere, so that 
the neutrino flux decrease with distance only for a dilution spatial factor, 

' Ru 



n^t - n^i = {n^t - ^fJq [—) (42) 

(for r > Ru) where {rii^^ — n-y^^ is the neutrino net number density at neutrinosphere, 
which we now pass to evaluate. Let us first concern on non electron neutrinos, that are 
produced entirely by thermal processes. The difference between neutrino and antineutrino 
distribution arises in this case only from the (slightly) different neutrino-nucleon elastic 
scattering cross section, and vanishes in the infinite nucleon mass approximation. The 
net number is then proportional to T /M^ {M^ is the nucleon mass) and is calculated in 



|T^; at the neutrinosphere we have 



{n^ - ^F )n - — ^ M^ T: (43) 

{x = fi, t) where S is a parameter depending on nucleon weak couplings, equal to 3.32 for 
neutrons and 2.71 for protons. 

The difference in z/g and Ve distributions induced by the different neutral current scattering 
is the same as in (|43|) ; however, for electron neutrinos and antineutrinos there is a more 
important contribution coming from charged current (elastic and inelastic) scattering. 
This induces a (not large) effective degeneracy parameter 77^ = fiu/T^, and the net electron 
neutrino number is approximately proportional to this parameter: 

(n,, - n^Jo ^ -V.TJ^ ■ (44) 

10 



Following |I|, Tjiy can be estimated from the fact that in a Supernova the ratio of the 
(average) energies in z/g and z/g +T'e, due to the deleptonization electron neutrinos from 
the core, is approximately equal to 5/8. From this we have 

r]^ ~ 0.29 (45) 

(but this result depends somewhat on T^, and on the average neutrino chemical potential 
in the core). 

From the above arguments we are now able to evaluate quantitatively the ponderomotive 
force exerted by neutrinos; we get 



-r-e _ _ -r-e^ ^ 
-6 



2.94xl0-«<! (^^l f^^^V - 3.2x10-2 (^^ f^i^ 
0.29/ K^MeV) 13.32/ V6Me\/ 



at) \t) ■■ —' <«) 



^\ f^^)' + 0.79 [l-\ f^ 
0.29/ V4Me\// 13.32/ V6Me\/ 



J^,? ~ - 1.5x10-^ <; l::^) I . ,:%. I +0.79 



Note that, for a given distance, the force acting on electrons is about a factor 13 greater 
than the one acting on protons, while electrons and neutrons experience about the same 
force. Observe also that while the contribution of v^ , v^ is subdominant with respect to 
that of z/g in T^ , each flavour of neutrinos contribute with about the same strength in 

'' V 1 "^ V ■ 

We can now estimate the maximum energy gain of the Supernova material behind the 
shock front due to the considered collective effects by assuming that electrons are rela- 



tivistic while the mean nucleon velocity is given by ~ JZT jM^q ^ 0.06. For K = 150 Km 



'^For illustrative purposes, we take the value 3.32 for S and consider e-neutrinosphere and x- 
neutrinosphere as located at about the same radius i?^, although with different temperature. This 
last approximation is justified from the fact that we will concern with effects induced by the ponderomo- 
tive force taking place not very near the neutrinospheres. Obviously, both these assumptions have to be 
relaxed in a detailed Supernova numerical simulation. 
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we clearly find that the electron energy gain is about a factor 3 x 10^^^ less than that of the 
standard results, and the electron heating is the major effect for collective interactions, 
contrary to the standard picture in which the nucleon heating (by inelastic scattering) is 
the most remarkable one. Moreover, in the proposed scenario, another suppression factor 
comes from the dependence on the angle between electron (or nucleon) motion and the 
pondero motive force (see Eq. (|4ll) ). In fact, the previous estimate of the heating rate 
assumes that nearly all the electrons (or nucleons) are polarized; then, although the mag- 
netic field contribution to Tv is unexisting (as established here), it would be nevertheless 
important for reaching the maximum effect 0, which is, however, completely subdominant 
with respect to the standard heating mechanism, the heating rate being some eV s~^ for 
r = Ry and decreasing as r~^ . 

Nonetheless, the action of the ponderomotive force due to neutrinos is not, in general, 
restricted to heat the material behind the shock (at least in a direct way). It could give a 
non negligible contribution to the theory of convection and instabilities during a Super- 
nova explosion. In fact we remind, from (^6|) -(48), that while Tl and T"^ are positive. 



Tl and JF" are negative. This means that e~ and "p are pushed in the outgoing radial 
direction, while e^ and n in the ingoing one, and this could be a relevant source of con- 
vection, as pointed out in [^, especially near the radiation bubble, where Ray leigh- Taylor 
instabilities already act. 

In conclusion, we have studied collective interactions of a beam of neutrinos/antineutrinos 
traversing a dense plasma of e^,p, n and applied the results to the case of a Supernova. 
We have found that the ponderomotive force exerted by neutrinos on the material behind 
the shock cannot substantially heat that and gives a negligible contribution to the revival 
of the shock for a successful Supernova explosion. In particular, the magnetic field present 
in the considered plasma does not contribute to the expression of the ponderomotive force, 
although it is relevant for the polarization effects induced by it in the medium. 
The collective interactions studied here can, nevertheless, play a relevant role in dense 



^In ||T2[ the mean electron polarization in a Supernova for the Wilson & Mayle model has been 
calculated. The authors find that for r ~ 100 -;- 200 Km it can be very high (near the unit value) allowing 
the magnetic field in this region to be as high as ~ 10^® Gauss. 
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stellar plasmas |T^ where they can provide an additional plasma cooling process through 
neutrino Landau damping of electron plasma waves, thus influencing the evolution of 
massive stars. 
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